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Abstract 

We consider the QCD corrections to the electric dipole moment of the neutron in the 
Minimal Supersymmetric Standard Model. We provide a master formula for the Wilson 
coefficients at the low energy scale including for the first time the mixing between the 
electric and chromoelectric operators and correcting widely used previous LO estimates. 
We show that, because of the mixing between the electric and chromoelectric operators, 
the neutralino contribution is always strongly suppressed. We find that, in general, the 
effect of the QCD corrections is to reduce the amount of CP violation generated at the 
high scale. We discuss the perturbative uncertainties of the LO computation, which 
are particularly large for the gluino-mediated contribution. This motivates our Next- 
to-Leading order analysis. We compute for the first time the order as corrections to the 
Wilson coefficients for the gluino contributions, and recompute the two- loop anomalous 
dimension for the dipole operators. We show that the large LO uncertainty disappears 
once NLO corrections are taken into account. 



1 Introduction 



CP violation plays a twofold role in SUSY model building. On the one hand, it is one of the main 
motivations to invoke New Physics (NP), since within the Standard Model (SM) it is not possible 
to construct a successful theory of baryogenesis, and also in the Minimal Supersymmetric Standard 
Model (MSSM) electroweak baryogenesis calls for additional sources of CP violation beyond the 
SM single phase in the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix [1]. On the 
other hand, CP violating processes provide very stringent constraints on NP. Indeed, the recent 
experimental progress in the study of FCNC processes allows us to conclude that most probably 
NP cannot contribute substantially (i.e. more than ~ 20 — 30%) to flavour and CP violation in 
AS = 2 and AB = 2 processes [2]. Since the MSSM contains tens of new sources of flavour and 
CP violation, this experimental observation is quite puzzling, leading to the so-called SUSY flavour 
problem, which is one of the main open issues in SUSY model building [3]. 

The new sources of CP violation present in the MSSM can be divided in two groups: the first 
one contains new phases that appear in flavour conserving quantities, the second contains those new 
sources of CP violation that are also new sources of flavour violation. While the latter are strongly 
constrained by K and B physics, at least for those flavour-changing parameters that connect the 
first generation to the other two,^ the new sources of CP violation that are not directly connected 
to flavour violation have little impact on FCNC processes [4] and are mainly constrained by the 
Electric Dipole Moments (EDM) of the electron and of the neutron [5-7]. A careful analysis of 
EDM processes is therefore mandatory in order to assess the allowed size of these NP contributions 
to CP violation and their possible effects in electroweak baryogenesis and in other CP violating 
processes. 

Surprisingly enough, while the study of FCNC processes in the MSSM has recently witnessed 
considerable theoretical advances, with the inclusion of Next-to-Leading Order (NLO) QCD correc- 
tions [8] and with the computation of the relevant hadronic matrix elements with Lattice QCD [9], 
not only no corresponding effort has been made in the study of the EDM of the neutron, but even 
incorrect LO results have been widely used in the literature. 

The present work aims to be a first step towards bringing the EDM analysis in SUSY at the 
same level of accuracy as the other FCNC and CP violation studies. In particular, we focus on the 
perturbative aspects of the QCD corrections to the neutron EDM. 

The paper is organized as follows. In section[21we discuss the neutron EDM at the LO in QCD: 
we give a complete LO formula for the EDM, correcting some errors in previous analyses, and 
we discuss the interplay of the various contributions. We then study the uncertainties related to 
the LO approximation, and we find that they are particularly large for gluino contributions. This 
motivates us to upgrade the analysis to NLO that is presented in the next section. We introduce 
the NLO QCD evolution and the NLO matching conditions for the gluino contribution. We show 
that after the inclusion of the NLO contributions the scale uncertainty is reduced down to ~ 2.5%. 
Finally we present some conclusions. 

^To get a feeling for numbers, the imaginary parts of squark mass terms connecting the first two generations are 
constrained to be < 10~^ — 10~* of the average squark mass. 
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2 Leading Order Analysis 



In this Section, we provide the full LO expressions for the neutron EDM in the MSSM, we dis- 
cuss the interplay of the various SUSY contributions and we study the uncertainties of the LO 
approximation . 



2.1 Anomalous dimension 

We write the relevant CP-violating effective low-energy Hamiltonian as 

Hcpv = C'MO'M + E + C^3(/u)03(m) , 



(1) 



where 



ol 



ol 



O3 = -Igsr^'Gl^GlPGl^e^^^^' 



(2) 



The index q runs over light quarks, and Qq = (2/3,-1/3) for up- and down-type quarks respectively. 
With this choice, all the operators have dimension six. 

Defining C = {Gl^C^^^C^) we write the renormalization group equation for the Wilson coeffi- 
cients as: 



dGjii) 
d Infj, 



l^G{^,) 



(3) 



where at the LO 7(0^) = (as/47r) 7'^''). 

Let us now discuss the LO anomalous dimensions of the operators in eq. ©. The anomalous 
dimensions of operators 0\ and Ol can be easily gleaned from that of the operators O7 and Og 
relevant in the b ^ s'j process (see ref. [10]). The anomalous dimension of the Weinberg operator 
O3 [11] and of its mixing with O2 was derived in ref. [12]. Therefore we get the following LO 
anomalous dimension matrix: 

/ 7e \ / 8C7^ \ 



7(0) = 



Iqe 1q 







8Gf WCp-'iN 







(4) 



V 7Gg 7g / V -27V N + 2nf + (3oJ 

where Gp = 4/3, N = 3, (3o = ^{UN — 2nf) with nj the number of active flavours. The Wilson 
coefficients at the hadronic scale can be easily obtained from those at a high scale from 

GfM = rj'^^Gf{f,s) + ^^{v''^-V^') Cliffs) + 

" 7g 



7e 



+ 



laqlqe V 



+ 



iGqlqe V 



GsM 
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C3(/"s) 
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(6) 
(7) 
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Table 1: Magic numbers Xf'- for the evolution from six to four flavours. See the text for details. 



where r] = asiiis)/as{fJ'H) and Ki = 7i/(2/3o). 

The operator basis in eq. is very suitable to discuss the anomalous dimension matrix. 
However, in order to avoid the explicit appearance of the strong coupling at the low scale in the 
operators, it is more convenient to introduce a slightly different operator basis 

Oa = -\r''G1pGiPGU^^^" (8) 

that defines our electric dipole (Oe), chromoelectric dipole (Oc) and Weinberg operator (Og) and 
whose corresponding Wilson coefficients can be easily obtained from eqs. by redefining the 

coefficients as follows: 

gsMCsifis) = V'^'^^^Cail^s). (9) 

To illustrate in a simple way the relation between the Wilson coefficients at the /xs scale and those 
at the fiH scale we take fis ~ nT-t and assume five fiavours of light quarks between the scales ^us 
and fiH, obtaining 

C!M = .ic|(;.s)+8(,i-,M)^+|l[l7,i-15,M-2,il^, (10) 

5 Q / 5 39 \ 

c!M = v^c!{t,s) + -(v^^-v^)Cait-is), (11) 
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CgM = V^Caifis). (12) 

It is interesting to note that in eqs. (|10H12p all the r/'s are raised to a positive power and then act 
as suppression factors. 

In general, SUSY masses are expected to be above mt while the hadronic matrix element is 
evaluated at a scale of the order of the neutron mass. In this situation it is more appropriate to 
consider the evolution from /is > mt to f^n < nT.^, i.e. from the six- to the four-flavour theory, that 
can be summarized via the so-called "magic numbers". In this case, the low-energy coefficients 
C{hh) = {Ci, Cc, Cg) are given in terms of the high energy ones as 

= E E X-^aMW''9s{f^sY'''''^'~'^Cj{^,s) , (13) 

j=l a,b=l 
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with Ya and Zj, given by: 




8 64 72 163 177 
75' 525 '175 '175 '175 



50' 25' 25 50' 25 



3 14 16 33 29 



} 



(14) 



and the nonvanishing entries in Xfj' are listed in Tabled These magic numbers have been obtained 
using the average values mt{mt) = 168.5 GeV, mb{mh) = 4.28 GeV and as{Mz) = 0.119. 

A comparison with previous evaluations of the LO anomalous dimension matrix is now in order. 
Several partial LO results are present in the literature [12-15] although the work of ref. [15], to be 
called ALN, can be regarded as the standard reference for the QCD correction to the neutron EDM 
with its numerical estimates of the QCD correction factors that have been and are still widely used. 
With respect to ALN our analysis differs in two aspects: i) we have included the mixing between 
the operators Oe and Oc that is neglected in ALN. ii) Our definition of the operator basis (eq. 
is different from that employed in ALN. In particular, we write explicitly in the definition of the 
operators Oe and Oc the mass of the quark, as well as in Oe the charge of quark, while in the 
operator basis of ALN the quark mass and charge is not present. Correspondingly the anomalous 
dimension matrix of ALN should differ from ours by a factor 7^ = —dCp. Taking into account this 
difference we find agreement with ALN in the anomalous dimension result for the chromoelectric 
and Weinberg operators. Instead, for the electric dipole operator we find that, with the conventions 
employed by ALN, the anomalous dimension should read 7e = —8/3 , i.e. it has the opposite sign 
with respect to the one quoted in ref. [15]. As a consequence, the QCD renormalization factor of 
the dipole operator, 77^^, that is estimated in ALN to be r/^^ = 1.53, should be rj^^ < 1 , i.e it 
does not enhance the CP violating effect but actually suppresses it. Employing the same values for 
strong coupling at the high and low scale used in ALN we get r]^^ ~ 0.61. 

In our view the definition of the operator basis we employ (eq. (jH}) has the advantage to make 
more transparent the behavior of the perturbative QCD corrections to the neutron EDM that in 
general give correction factors that decrease the amount of CP violation generated at the high scale. 
In the ALN operator basis this effect is somewhat hidden by the fact that the quark mass entering 
their Wilson coefficients has to be taken at the high scale and mq{ps) < n^qifJ-ii)- It should be 
noticed that the dependence of the Wilson coefficients upon the quark mass can also appear in an 
indirect way, e.g. through the matrices that diagonalize squark masses. 

2.2 Hadronic Matrix Elements 

In order to compute the EDM of the neutron the matrix elements of Oi, Oc, Oq between neutron 
states are also needed. At the moment a result from Lattice QCD is not yet available, although 
first steps in this direction have been recently made [16]. Several alternative approaches have been 
used to estimate these matrix elements, as QCD sum rules [17] or chiral Lagrangians [18]. In this 
paper we are mainly concerned about perturbative aspects of the EDM calculation, thus we are 
going to use the simplest estimates of the operator matrix elements. In particular, for the electric 
dipole operator, we use the chiral quark model where the neutron is seen as a collection of three 
valence quarks described by an SU(Q) symmetric spin-flavour wave function. In this model the 
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neutron EDM is related to that of the valence quarks by 



(15) 



where 



dq = eQgmg{fiH) C| {fin 



) 



(16) 



is the quark EDM. Concerning the contribution of the chromoelectric and Weinberg operators to 
dn the simplest estimate is obtained via naive dimensional analysis [19] giving 



where A ~ 1.19 GeV is the chiral-symmetry-breaking scale. We notice that in eqs. ()16() and H17|) 
ruq is computed at the hadronic scale, while in the expressions for the Wilson coefficients at the /is 
scale the masses, as well as gs, are computed at the high scale. 



The discussion in the previous sections has been general and applies to any model in which CP 
violating effects are generated at some high scale. In this section we focus on the minimal su- 
persymmetric standard model (MSSM) with complex parameters assuming also that the trilinear 
SUSY-breaking scalar couplings are proportional to the corresponding Yukawa coupling. In this 
model besides the SM CP violating phases {5c km, (^qcd) that will be neglected in the present 
discussion as well as the flavour mixing, there are new phases associated to the /x-term in the 
superpotential, the supersymmetry-breaking parameters of the gaugino mass, the trilinear scalar 
couplings A, and the bilinear scalar term in the Higgs potential. However not all the phases of these 
quantities are physical. It is then possible to assume the gaugino masses and the bilinear term to 
be real so that we are left only with two CP violating phases, one associated with the fj, term, 
{4>fj,), and the other with the trilinear scalar coupling, {(Pa), that is in general flavour-dependent. 
These phases will be present in the mass matrices of squarks, charginos and neutralinos inducing 
an EDM at the quark level. In particular in the squark mass matrix 



the only complex parameter is the left-right mixing term Xu = ^* — ficot (3, = A'^ — /itan/3, 
where both phases are present. In Xq, tan/3 = V2/V1 is the ratio of the VEVs of the Higgs fields. 
Instead in the chargino mass matrix 




(17) 



(18) 



2.3 Wilson Coefficients in the MSSM 




(19) 




(20) 



as well as in the neutralino one 




/ Ml 

M2 

-C/3 Sw Mz C(3 cw Mz 

\ 3/3 SW Mz -8/3 CW Mz 



-Cf3 Sw Mz Sf3 sw Mz \ 
Cf3 Cw Mz -S/3 Cw Mz 

-/i 



(21) 
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the only complex term is ^ and therefore only cp^ is present. In eqs. (|2UI21|) sp = sin/3, cp = cos/3, 
sw, cw are the sine and cosine of the weak mixing angle and Mi^2 are the soft-breaking gaugino 
masses associated to the U{1) and SU{2) groups. 

To present our results in a transparent way we perform the computation of the Wilson coef- 
ficients using current eigenstates for squark fields. In this basis the squark propagator is a 2 x 2 
matrix given by 



(fc^ — m?^) (/c^ — 77i?^) y mqXq k"^ — m?- 



(22) 



92 



where rrtq^ are the eigenvalues of M|. Neglecting Oim^) terms we have that rriq^ ~ f^qL-> 
so that Ag(A:) reduces to 



A,-(fc) 



mqX* 



(23) 



k'^—m- 
12 



We notice that within this approximation the left-right propagator is still exact. 

The Wilson coefficients of the operators Oe, Oc are generated at the one- loop order while that 
of 0(3 appears for the first time at the two- loop level. At the LO we can then set Cdf^s) = and 
write 



ex 



ex" 



(24) 



and similarly for Cc- 

We find for the gluino contribution 



47r m? 

Qs a-s 



Im 



m-gj 3 



X, 



Alii — C (xi,X2) 
47r m| \mh 



(25) 
(26) 



where Xi = m^/rn^_. The explicit expressions for the functions B and C as well as those entering 
in the chargino and neutralino contributions are collected in the Appendix. 
For the chargino we have 



Q-(/^s) 



C^4-(/^s) 



1 2 



E 



1 



— Im 

1 2 



m 



Mw 



{lm(y, 



i2 Uil 



{Qd-Qu)A{yi,) + QdB{yl 



2 

\ m _ 

\ Xi 



Vi2 Ui2 Yd 



{Qd - Qu) hy\ , 4 ) + Qd B[y\ , ) 



(27) 



m 



Mw 



{Im {Vij U,2) [{Qu - Qd) A{yl^) + Q„ i3(4)] 



C^cV(/^s) 



Im 



V^^2 Ui2 



9s a 



E 



^ V2S/3 



iQu-Qd)A{yi^,ylJ + QuB{yl,yi 
{lm{V,2Ua)B{y'^J 



(28) 
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Vi2 Ui2 



C4-(/^s) 




(29) 




(30) 



where y 



m 



?_/m?_^. , U and y are the matrices that diagonahze according to U*M^,V ^ 



and d are the Yukawa coupUngs of the up and down quarks in units of e/ sw- In eqs. l(77HSn|) 
we have also written expUcitly the contributions proportional to the mass and the Yukawa coupling 
of the light quarks to show that the phase combination that enters in the gluino contribution is 
actually present in the chargino term in a suppressed way. Indeed, in the chargino contribution 
the only relevant phase is cj)^, hidden inside the matrices U and V . This can be seen explicitly in 
the following simplified expression, obtained by neglecting the contributions proportional to quark 
masses and Yukawa couplings: 



1 



Im 



Mw 



Analogous expressions can be obtained from eqs. 

Im(V-2C/a) 



Im {Vii Ui2) 



{Qd-Qu)A{yl)+QdB{yl 



EHI)-(|Sni) with the substitutions: 
> Im(C/^C/a/i) , 



(31) 



(32) 



Finally the neutralino contribution, neglecting terms proportional to the the quark masses, is 
given by: 



a 



E 



1 



2 

Xi 



Im 



+ 



Im {K^,K'^,) 



9s C'g^o 



(33) 
(34) 



where = m-n/m? , and 
3 xr ij 



Ul 



V2 



1 



V2 



z 



ii 



1 



\/2tan6lH/ Qd Z'^ 
1 



V2cp 



Z 



i3 



Qrf + - ) tan 9w Z 



il 



2 



(35) 

(36) 
(37) 

(38) 

(39) 
(40) 
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As can be 



In eqs. (|35I40|) Z is the matrix that diagonalizes according to Z*M^Z 
seen from eq. (jHSJ in the neutrahno contribution both (/>^, through the matrix Z, and the Xq phase 
combination are actually present. 

We noticed that the results reported in eqs. ((^ - (|5n|) and (|5S)) - H4Up are fully in agreement with 
those in ref. [5,6] and represent the lowest order approximation. 



2.4 LO results 

In this section we investigate, at the LO, the effect of QCD corrections on the neutron EDM, 
to assess whether they can significantly reduce the individual gluino, chargino and neutralino 
contributions making the EDM constraint on SUSY phases less severe. 

To discuss in a simple way the effect of the QCD corrections and in particular the importance 
of the mixing between Og and Oc that was neglected in previous analyses, we consider eqs. (|1UI12() 
assuming = 0.3 and setting Cq = 0. Then 



0.43C|(^s) 

0.88C,^(^s) 
0. 



0.38 



Cliffs) 



(41) 

(42) 
(43) 



Thus, if Ci{fis) — C'c (ms)/5s(/"s) the resulting Cid^n) is strongly suppressed. With our definition 
of the operators the above situation is achieved when the gluon in a diagram contributing to Cc is 
attached to a squark of the same charge of the external quark. This case is realized in the neutralino 
contribution (see eq. (EU)- Indeed, we can estimate the neutralino contribution to d^, d^, as given 
in eqs. (|15I17() . by employing the Wilson coefficient at the low scale evaluated via eqs. P0I12() with 
5s(/is) = 1.22. We get 



^rudM ( -^0.05 ) C^oifis) - muM ( ^0.05 ) &(/is) 



0.02 C^^o (ms) - 0.01 C^A^'s 



ex" 



0.88— 
e0.08 



rriuiliH) C"^o(/is) + mdiiiu) C^^fis] 



ex" 



(44) 



(45) 



Eqs. H44I45|) show that the individual quark EDMs are strongly suppressed by the QCD corrections. 
A so large effect is actually specific to the neutralino contribution because of the simple relation 
between C^-o and C^-o (eq. (jBU)- The general case is more complicated and the resulting effect 
depends upon the relative sign between C|(/is) and Cc(/is)- 

It is not our purpose here to perform a general analysis of EDM constraints on SUSY models. 
Rather, our aim is to illustrate the impact of QCD corrections on the computation of the EDM. To 
do so, we study a specific point in the SUSY parameter space that we choose with a mass spectrum 
similar to that of the benchmark point la of the Snowmass Points and Slopes ^ (SPS) [20]. The 
SPS benchmark points are actually defined assuming real parameters, however we take the mass 

^We use the low-energy spectrum obtained from the tool [21] available on the Web at the URL: 
http : / /kr ami . home . cern . ch/kraml/ comparison/ 
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Figure 1: Gluino contribution to the neutron EDM versus (p^ (top) and (pAu = <pAd (bottom), 
without (left) and with (right) QCD corrections. The solid line is the d^ contribution, while the 
dotted one the corresponding contribution. 

spectrum of point la as indicative, also in the case of complex parameters, of possible mass values 
of a "typical" mSUGRA scenario with an intermediate value of tan/5. We take rUg = 585 GeV, 
Ml = 100 GeV, M2 = 190 GeV, rui,^ = 540 GeV, m^^ = 525 GeV, mj^ = 550 GeV, mj^ = 520 
GeV, 1^1 = 355 GeV, tan/3 = 10, \Ad\ = 855 GeV, = 675 GeV and /is = 465 GeV. To obtain 
dn we have chosen the hadronic scale fin = 2 GeV with m^{fiH) = 7 MeV and rUuif^H) = 3 MeV. 

In fig. ^we show the effect of QCD corrections on the gluino contribution to the EDM of the 
neutron. In the figure we plot the absolute value of dn as a function of (j)A and (/>^ with the other 
SUSY parameters set to the values listed above. As can be seen comparing the plots on the left, 
which are obtained without QCD corrections, to the ones on the right, where QCD corrections 
are included, the effect in this case amounts to ~ 10%. Fig. Hand fig. 01 show the corresponding 
analysis for chargino and neutralino contributions, respectively. Fig. [21 shows that for chargino 
contributions the inclusion of QCD corrections reduces the amount of CP violation generated at 
the Us scale by a factor ~ 50%. Finally, the simple analysis of the neutralino contribution discussed 
above is substantiated by fig. O where this strong reduction is clearly visible. 



10 




Figure 2: Chargino contribution to the neutron EDM versus (pfj,, without (left) and with (right) 
QCD corrections. We show only the contribution because the one is negligible. 




Figure 3: Neutralino contribution to the neutron EDM versus (j)^ (top) and 4>Au — 4'Ad (bottom), 
without (left) or with (right) QCD corrections. The solid line is the contribution, while the 
dotted one the corresponding d^ contribution. 
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A popular mechanism [5] invoked to suppress the neutron EDM without resorting to extremely 
small phases or very heavy SUSY particles is the search for regions of the parameter space where 
cancellations among the three different contributions are active. It is always possible to find regions 
of the parameter space where contributions depending upon different parameters cancel each other, 
although it can be questioned if these regions can be representative of general situations. With 
respect to this, it is interesting to note that, since the neutralino contribution is always much more 
suppressed by QCD corrections than the gluino and chargino ones, the cancellation mechanism 
among different contributions invoked in ref. [5] should actually work between the gluino and 
chargino only. However, these two contributions depend upon different phase combinations. As an 
example, (pA is only present in the gluino contribution. 

2.5 Uncertainties of the LO analysis 

In the above analysis all the uncertainties of the LO computation have been neglected. The uncer- 
tainties connected to the nonperturbative evaluation of the hadronic matrix elements go beyond 
the scope of this work, since we focus our analysis on the perturbative aspects of QCD effects. 
Therefore, let us assume that some nonperturbative method such as Lattice QCD will produce in 
the future the necessary matrix elements at a scale = 2 GeV, so that we fix the hadronic scale 
in our analysis. Then, we are left with the uncertainties connected to the matching between the 
full and the effective theory at the scale /is- 

It is well known that in the RGE improved perturbation theory there remain unphysical //s- 
dependences which are of the order of the neglected higher order terms. Usually, this uncertainty 
can be estimated by varying the matching scale in a (arbitrarily chosen) given range. However, for 
the EDM computation, there are further sources of uncertainty. All contributions depend upon the 
squark masses, but the precise definition of these masses cannot be fixed at LO, so that one can 
use pole, DR or any other squark mass. Indeed, the difference between the results obtained using 
two different mass definitions is of higher order in and provides an estimate of this additional 
LO uncertainty. Furthermore, the gluino contribution also depends on the gluino mass and, more 
important, on the strong coupling. Neither the definition of the gluino mass nor, in principle, 
the scale of is fixed at LO, so that they constitute another source of uncertainty. All these 
uncertainties can be ameliorated only by a NLO calculation. 

In fig. 0] we illustrate the LO uncertainty only due to the choice of the matching point. Here 
and in the following NLO analysis we will use an average squark mass m? = (m?^ + m'^^)/2. In the 
figure we plot the LO gluino (right) and chargino (left) contributions to \dn\ — as a function of 
the matching scale with and DR masses evaluated, for simplicity, at the scale /is- As expected 
the /is dependence is more pronounced in the gluino case and amounts to 10-15 % while in the 
chargino contribution it reaches at most 4 %. 

The LO gluino contribution shows a substantial uncertainty that, including all effects, can be 
expected to be ~ 20%. To reduce it to a level comparable to that of the chargino contribution one 
needs the NLO computation of this contribution that will be discussed in the next section. 
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3 Next-to-Leading Order Analysis 

In this section we do not attempt to perform a complete NLO analysis of the QCD corrections to the 
neutron EDM, instead we focus on the relevant pieces needed to discuss the reduction of the scale 
dependence of the gluino contribution. We present the NLO anomalous dimension matrix for the 
electric and chromoelectric operators and the NLO Wilson coefficients for the gluino contribution. 
For completeness we present also the Wilson coefficients of the Weinberg operator. We recall that 
at the LO Cc{^h) = 0, therefore to obtain the NLO result it is sufficient to know the LO 7g, 7g(j 
entries of the anomalous dimension matrix. 



3.1 NLO anomalous dimension 



The discussion of the NLO anomalous dimension matrix is more easily accomplished in the O^-Oa 
basis of eq.Q. Indeed in this basis the anomalous dimension matrix can be organized in powers of 
as as 



7 



4tt V47r 



7 



(1) 



where 7^*^^ is given in eq.(|lj) and 



/ (^^N - 16 Cf 



'^N-32Cf 



fuf) Cf 



fuf) Cf 



458 
9 



^ + — 





214^2 _|_ 56_IV 



13 



9 TV 



Nnf 



(46) 



(47) 



^ * * * y 

We have computed the NLO anomalous dimension in ea. ()47() and our results are in complete 
agreement with those obtained from the known NLO anomalous dimension of the O7 and Og 
operators in the 6 — > 57 process [22]. 

The entries in the third row are unknown but not needed at the NLO, since they contribute only 
to the (z3) sector of the NLO magic numbers, and the initial coefficient of the Weinberg operator 
is vanishing at the LO. The corresponding Wilson coefficients in the O^-Og basis can be easily 
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ab 


43 


43 


54 


55 


63 


ij 


11 


12 


12 


22 


12 


'J 


11.301 


85.158 


-79.353 


9.9191 







-8.7762 


-70.209 


65.693 


-6.9887 







-0.39785 


-3.1828 


3.7203 


-0.65874 


-0.46504 



Table 2: Magic numbers R1° , W^f and SX^^^ for the NLO evolution from six to four flavours. 



obtained using eq. @ . The simplest way to present the NLO evolution of the Wilson coefficients is 
via the magic numbers. Referring to the O^-Og basis we can write for a generic scale // 



c(M) = c(°)(/.) + ^c«(^) 



(48) 



where cf^\iJ,s) is the LO Wilson coefficient at the scale and its evolution from > i^t to 
IJ,H < rub is given by ea. ((T^ . The evolution of C^^\fi) from /is > mt to fin < can be 
summarized in the following way 



j=l a=l 6=1 



(49) 



The relevant entries of the NLO magic numbers SXf^ , Wj" and i??- are given in Table |2 As 
expected, the evolution of C^^^(/xs) is dictated by the magic numbers derived in the LO case (see 
eq.®). 



3.2 NLO Wilson coefficients 

The computation of the matching conditions at the NLO level can be divided in two parts: the 
matching conditions for the helicity flip operators, Oe, Oc, and those for the Weinberg operator 
Oq. Concerning the latter, the light quarks give a vanishing contribution so that the only relevant 
contribution is due the top quark that has no associated effective theory to subtract and therefore 
no infrared (IR) divergent terms to deal with. Instead, for helicity flip operators the computation 
of the NLO matching condition is, in general, a very complicated task. At the moment, even for 
a process like h 57, that has been investigated in great detail in the last ten years, we have 
not yet obtained in the MSSM the complete NLO matching conditions but only partial results are 
available [23,24]. 

We begin discussing the NLO matching conditions for Og, Oc- In the actual computation two 
strategies are at hand. One can match matrix elements of operators belonging to a basis, like the 
one in eq.(jSJ, obtained enforcing the equation of motion, a procedure that however requires, in 
general, an asymptotic expansion of the relevant diagrams in the external momenta. Alternatively, 
one can use a larger off-shell basis and perform the matching on the off-shell matrix elements. In this 
case, one can use the freedom of the off-shell status to choose a suitable kinematical configuration 
such that the relevant Feynman diagrams can be evaluated using ordinary Taylor expansions in the 
external momenta. The latter strategy, applied in ref. [25] to the NLO matching conditions of the 
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Ai3 = 1 magnetic and chromomagnetic operators, has been employed by us in the present work 
fohowing closely ref. [25] to which we refer for technical details. 

The off-shell operator basis relevant for our calculation is obtained by supplementing the basis 
in eq.® with the two operators [26] 

= -\qjpV>im (50) 

where is the SU{3) x U{1)q covariant derivative. The relevant terms in the Wilson coefficients 
are extracted via the use of the projector 

+ ml [(n-l)(4 ^^i(2pi -(7)'^-4g2pA') + 2n(72g^_2gV ^] } 75 (51) 

assuming an off-shell kinematical configuration defined p\ = = Pi ' Q = l'^ where pi 
and q are the momentum carried by the incoming quark and the external boson, respectively, 
and <C m^. The projector works by contracting it with the amplitudes and taking the trace, 
so that in eq. (|51|) fi is the index carried by the external boson while n is the dimension of the 
space-time. The use of an off-shell kinematical configuration induces in the result for the "full" 
theory the appearance of terms that behave like I/jti^, Inm^ as mq 0. These infrared terms are 
eliminated by corresponding terms in the effective theory once the off-shell basis (Of, Oc, Om, 0'^2) 
is employed so that the Wilson coefficients are free of infrared terms as they should. 

To simplify the calculation we compute the NLO gluino contribution to the matching conditions 
retaining only one source of CP violation, namely we keep only one power of Xq, discharging terms 

with n > 1. We also work in the limit of = rriq^ = rriq with rriq common to all squark 
flavours and taking all quarks massless but the top one. Within this framework the NLO gluino 
corrections can be written as 



clf\m) 




Fi (xg) + 4F2 (xg) + F2 {xg, xt) + Re ( I {xg, xt) 



m- 



+ Im(^l Re(^] N2ixg,xt) \, (52) 



^ 47rmg yiTT-gJ 



Fs {xg) + 4F4 (xg) + F4 {xg,xt) + Re ( I {xg, xt) 



+ Im(^!^jRe(^)iV,(,-,,,)| (53) 

where in the above equations the upper line represents the CP violation induced by the left-right 
entry in the mass matrix of the squark of type q while the lower one the corresponding effect due to 
the stops. We have further divide the former contribution into the the part due to the quark and 
squark of type q, that of the other four squarks and massless quarks (the first two terms), and that 
due to the top and stops including the mixing (the last two terms). In eqs. (|52l53|) Xg = m^/m?-, 
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Figure 5: The fig scale dependence of the gluino contribution to d^ at LO (upper curve) and NLO 
(lower curve). 



Xt = m1/m?^, , where the gluino and squark masses are assumed as DR parameters, and 

F, = Gi + Ai In , hm Fi{xg,xt) = Fi{xg). (54) 
The exphcit expressions of the functions Gi , Aj, Ni are reported in the Appendix. Defining 

Ai 



-^Im ) (Ai + 5A2) , 



Ao 



Qs as 



Im(^)(A3 + 5A4) 



we have that the coefficients of the ln(^g/m?) terms satisfy 



A. 



^ , (0) 



k=g,q 



druk 



(55) 



with 



SUSY 



3N-nf , 7W=2(3iV-n^) , j!^>= ACp ml/ml , 7^" 



(0) 



-g/ ■■''q 



„(o) 



-2Cf 



guaranteeing the cancellation of the dependence to 0{a1) in eq.(@SJ. We observe that the effect 
of the term in the square brackets in ea. H55() is to shift the coupling and the mass parameters 
appearing in cf'^ from the scale /xs to niq. 

As we anticipated in sec. 12.51 the inclusion of the NLO matching for the gluino contributions 
reduces the perturbative uncertainties down to a completely negligible level. Moreover the inclusion 
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of NLO corrections produces a non-negligible effect. In fig. 1^1 we plot the scale dependence of the 
gluino contribution at the LO (upper line) and at the NLO (lower line) level. As shown in the 
figure the inclusion of the NLO contribution greatly reduces the scale dependence of the gluino 
contribution and lowers \dn\ of about 10%. 

Finally, for completeness, we consider also the Weinberg operator. At the scale two-loop 
diagrams where top and stops together with the gluino are exchanged contribute to C^'^ . The 
relevant expressions can be gleaned from ref. [27] obtaining 

4i^(Ms) = f^Imf^) H{x,,x,), (56) 

where the function H is found in the Appendix. However, when the evolution down to a four-flavour 
theory is considered one has to take into account also the shift in Cq'^ induced by the Oc operator 
at the mf, threshold [12, 13, 27] or 



where 



= ^^C'glMs), (58) 
C'}^\mb) = Vti^Clfi^s) (59) 



with rib = as(/is)/as("^fe)• 



4 Conclusions 

In this paper we have discussed the LO and NLO QCD corrections to the electric dipole moment 
of the neutron in the MSSM. We pointed out the importance of the mixing between the electric 
and chromoelectric operators that was always neglected in previous analyses. Also we noticed that 
the QCD renormalization factor of the dipole operator in absence of mixing should be less than 1 
while its widely used estimate is tj^^ = 1.53 [15]. 

In the MSSM the prediction for the EDM can easily clash with the experimental upper bound 
\dn\ < 6.310-26 e-cm [28] if the phases in the mass parameters are arbitrarily chosen. To avoid 
conflict with the experimental bound one can consider models with approximate CP symmetries [29] 
or flavour off diagonal phases [30] where small phases can be naturally obtained. Another possibility 
is represented by a cancellation mechanism among different contributions to the quark EDM which 
could allow 0(1) phases. However we noticed that, because of the mixing between the electric and 
chromoelectric operators, the neutralino contribution is always much more suppressed than the 
gluino and chargino ones so that the cancellation mechanism should actually mainly work between 
the latter contributions that, however, depend upon phases connected with apparently unrelated 
terms in the MSSM Lagrangian. 

Our results show that the NLO corrections we considered lower the prediction of the EDM 
of about 10% with respect to LO. Moreover the dependence on the matching scale is drastically 
reduced. Clearly, a complete NLO analysis will require, besides the gluino contribution we focused 
on, the other two contributions, in particular the chargino one. 
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We have mainly focused on the perturbative aspects of the problem but, in order to achieve a 
complete analysis of the SUSY constraints from the neutron EDM, a lattice computation of the 
relevant matrix elements is mandatory. 

Finally we notice that our analytic formulae for the magic numbers can also be used for the 
evolutions of the Wilson coefficients in any extension of the SM with new heavy particles, unless 
large four fermion operators are generated at the matching scale. This happens, for example, in 
supersymmetric models with large tan/? [31]. 
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Appendix 



In this appendix we give the explicit expressions of the one- and two-loop functions that appear in 
the Wilson coefficients. 

The explicit expressions of the one- loop functions entering in eqs. H24II34() are given by 



A{xi,X2) 
B (xi,X2) 



xiX2 {xiX2 -I- a^i -I- a;2 - 3) 



XiX2 



2fxi 



1) 



[X2 



11 



{Xi - X2) 



xi In xi X2 In X2 



(Xi - 1)^ (X2 - 1 



X1X2 (3xiX2 — Xl — X2 — 1) X1X2 



2(xi 



1)' (^2 



11 



+ 



(xi - X2) 



' (x 
xf In xi X2 In X2 



Ixi 



1)' (^2 - 1) 



X1X2 (3xiX2 + 5x1 + 5x2 — 13) 



3(xi 
X1X2 



1)' (^2 



11 



A{x) 
B{x) 



3 (Xi - X2 

X (x — 3) 
2(x- 1)^ 
X (x -|- 1) 



xi (xi — 9) In xi X2 (x2 — 9) In X2 



+ 



X Inx 
x^ Inx 



(^2 - 1)^ 



2(x-i)^ {x-iy 

In the case of equal masses the functions A, C reduce to 



A{x) 
B{x) 
C{x) 



x^ (x -I- 5) x^ (2x -1-1) Inx 
~2 (x- 1)^ ^ (x - 1)^ 
x^ (5x -I- 1) x'^ (x -1-2) Inx 



2 (x- 1)^ 
2x2 {x + ll) 
3 (x - 1)^ 



+ 



(X - 1)^ 
x^ (x^ — 16x — 9) In 



X 



3(x-l)' 



(Al) 



(A2) 



We list now the two- loop functions that appear in eqs. (|52l53j) . To simplify the expressions, we 
perform an expansion in the top mass reporting only the first term in xt- In particular, recalling 
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eq. (|54|) . we write 



~ G2 {Xg) + Xt G\ {X~g) + A2 {Xg) 111 



~ G4 (xg) + G\ (xg) + Xt In Xt 5*4 (xg) + A4 (xg) 111 
iVi(xg,xt) ~ ./x'tN- (xg) + y/TtlnxtRiixg) . 



2 
3. 
2 

2 
3. 



We find 

8x2 (51x3 + 413x2 - 1473x - 251) 16 (Sx^ + 293x - 13) 

Gi X = ^ 2 ^ J ^Li2 1-x 

^ ^ 27(x-l)^ 27(x-l)^ ^ ^ 

8x2 (i27x'' - 1075x3 + 480x2 + 405x + 27) ln2 x 
^ 27(x- 1)*^ 

8x2 (^433.4 _^ 228x3 - 1105x2 - 548x + 8I) Inx 

27 (x- 1)^ 



N{ (x) 



(A3) 



64x2 (x2 + 4x - 2) 16x2 (x + 5)^ , , 16x3 (^2 + 4x - l) Inx ^^^^ 

G2{x) = ^ -4 3^ Li2(l-x) (A4) 

3 (x — 1) 3 (x — 1) (x — 1) 



X- (129x3 - 2903x2 + 1083x + 21851) 4 (113x2 + 281x + llO) 



x2 , , ^ , , 

G3 (x) = ^ 1 Li2 (1 — x) 

^ ^ 54(x-l)' 27(x-l)^ ^ ' 

^ x2 (539x^ - 2282x3 - 6744x2 + 7578x + 62l) ln2 x 



27{x-iy 

X- {96x^ - 5019x3 + 13357x2 + 10583x + 1719) Inx 
54 (x- 1)' 



+ ^ 77-. :T5 (A5) 



x2 (17x2 - 310x + 101) 2x2 (2; - 40) , , x2 (2x3 - 67x2 + 4x - 3) Inx 

G4 X = ^ J ^ + ^ ^ Li2 (1 - X + — ^ = '- 

' 3{x-lf 3(x-l)^ ^ ^ {x-lf 

_ 8x2 (^2 ^ ^ 175) 3^^2 + 3) 32x3 (2^ + 7)1^^ 

^^^''^ - ^(^^Ij^ + (x-1)^ Li2(l-x)- (A6) 

, , , x2 (17x2 + 230x - 4387) 2x2 2 _ g _^ 

Giix) = ^ 2 ^ + ^ , ^Li2 1-x 

' 18(x-l)' {x-lf ^ ' 

x2 (53x2 - 305x + 18)lnx , , 

+ — i s Inx A7) 

3(x-l)^ ^ ^ 

x2(x2 + 10x + l) 6x3(2; + l)lnx 
64 (x) - - — ^ \ - —5 (A8) 

(x - 1) (x - 1) 



4xt (x3 + 8x2 _^ 173^ _^ 34) 32^1 + 2) In : 



9(x-l)^ 3(x-l) 
16xt (5x2 _^ I'j^ _^ 2) 



3(x- 1)- 



Li2 (1 - x) (A9) 
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^2(^) = - :I - ;^ 74 Li2(i-x) (mo) 



Niix) 



4xt {x^ + X + 34) 16xi (x + 2) ^ 

9(x- 1)^ 3(2;-l)^ 
xi (llx^ - 158x2 + 1225x + 542) 2xi (2x2 - 86x - 9) Inx 
18(x- 1)^ ^ 3(x- 1)^ 

2xt (5x2 _ ;^g^^ _ 52) 

+ 4? T^s (1 - x) (All) 

3 (x — 1) 

xt (11x3 _ 114^^2 Y53x - 758) 6xi (6x - 1) In x 
18(x- 1)^ (x- 1)^ 



2x2 (x2 - 89x + 52) 
3(x - 1)' 



+ — 775 ^Li2 (1 - x) (A12) 



Rlix) = 8x-17)_3xi(3x + l)lnx ^^^3^ 



^ ^ 2 

X2 X 



2(x-l)^ (^-1)^ 
Rlix) = (^^-8x-17)_3xi(3x + l)lnx ^^^^^ 



2(x-l)^ (^-1)^ 

16x2 (8x3 + 13x2 - 176x - 37) 16x2 (29x3 - 97x2 - 115x - 9) Inx 
Ai X = ^ J ^ ^ . ^ A15 



, , 8x2 (7x2 + 16x + 1) 16x3 (x2 + 7x + 4) Inx 

A2(x) = 4- -^^ (A16) 

3 (x — 1) 3 (x — 1) 

8x2 (2x3 - 83x2 + 268x + 197) 4x2 (59x3 - 67x2 - 643x - 117) Inx 

= ^ ^7 — T^^ ^ ^7 — — 

9 (x — 1) 9 (x — 1) 



^4 (a;) = ^— -4 ^ + ^ -5 (A18) 



2x2 (x2 + 64x + 31) 2x2 (x3 - 29x2 - 59x - 9) In; 

3(x-l)^ 3(x-l)' 

where Li2(z) = — [ln(l — is the dilogarithm function. 

We observe that we have to take care of the fact that the naive dimensional regularization (NDR) 
we used violate super symmetry, because the gauge boson and gaugino interactions with matter are 
different at one loop. Supersymmetric Ward identity can be restored with an appropriate shift 
in the gluino-squark-quark coupling and with a shift of the gluino mass [23,32]. Explicitly, the 
coupling and the gluino mass in the lowest order formula (eqs. (|25|26|) ). must be replaced with 

OL 4 

m-g -> my{\ + ^3) . (A20) 

Finally, the function H entering in the Weinberg operator (eq. H56|)) is given by 

H{xg,xt) ~ H (xg) +xtHl{xg) + xtlnxtH2{xg) 

_ x2(x + ll) x2 (x2 - 16x - 9) Inx 

^^"^^ ~ 3(x-l)^ ^ 6(x- 1)^ 

X (5x3 + 265x2 + 455x + 27) 2x2 (4x3 - 93x2 - 258x - 8l) Inx 

(x) = — ^ = i- H t: 

^ 6(x-l)^ 3(x-l)^ 

2x2 _ i2x^ - 51x - 18) 
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X 



(llx^ - 223x2 - 259a; - 9) (x^ - 12x2 - 51x - 18) In; 



6(x- 1) 



+ 



(x-1) 



(A21) 
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